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Abstract 

We obtain Lewenstein-Sanpera decomposition of iso-concurrence decomposable states of 2(g) 2 
quantum systems. It is shown that in these systems average concurrence of the decomposition 
is equal to the concurrence of the state and also it is equal to the amount of violation of positive 
partial transpose criterion. It is also shown that the product states introduced by Wootters in 
[W. K. Wootters, Phys. Rev. Lett. 80 2245 (1998)] form the best separable approximation 
ensemble for these states. 
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1 Introduction 

During the past decade an increasing study has been made on entanglement, although it was 
discovered many decades ago by Einstein and Schrodinger |Q, [j. This is because of potential 
resource that entanglement provides for quantum communication and information processing H, ||, 
[|. Entanglement usually arises from quantum correlations between separated subsystems which 
can not be created by local actions on each subsystems. By definition, a bipartite mixed state p is 
said to be separable if it can be expressed as 

p = ^ Wi pf' pf ) , wt > 0, ^2wi = l, 

i i 

where pf and p\^ denote density matrices of subsystems 1 and 2 respectively. Otherwise the state 
is entangled. 

The central tasks of quantum information theory is to characterize and quantify entangled states. 
A first attempt in characterization of entangled states has been made by Peres and Horodecki family 
[||, |?|] . It was shown that a necessary condition for separability of a two partite system is that its 
partial transposition be positive. Horodeckis showed that this condition is sufficient for separability 
of composite systems only for dimensions 2 % 2 and 2 % 3. 

There is also an increasing attention in quantifying entanglement, particularly for mixed states 
of a bipartite system, and a number of measures have been proposed (5|, H ||, 10 1. Among them the 



entanglement of formation has more importance, since it intends to quantify the resources needed 
to create a given entangled state. 



An interesting description of entanglement is Lewenstein-Sanpera decomposition [11|. Lewen 



stein and Sanpera in [11] have shown that any two partite density matrix can be represented 



optimally as a sum of a separable state (called best separable approximation (BSA)) and an entan- 
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gled state. They have also shown that for 2-qubit systems the decomposition reduces to a mixture 
of a mixed separable state and an entangled pure state, thus all non-separability content of the 
state is concentrated in the pure entangled state. This leads to an unambiguous measure of entan- 
glement for any 2-qubit state as entanglement of pure state multiplied by the weight of pure part 
in the decomposition. 



In the Ref. [11], the numerical method for finding the BSA has been reported. Also in 2 ® 2 



systems some analytical results for special states were found in [12]. An analytical expression for 



L-S decomposition of Bell decomposable states and a class of states obtained from them via some 



LQCC action is also obtained in [13] 



In this paper we introduce iso-concurrence decomposable (ICD) states. We give an analytical 
expression for L-S decomposition of these states to show that the average concurrence of the de- 
composition is equal to their concurrence in these states. As a byproduct we also show that the 
concurrence of the considered states is equal to amount of violation of positive partial transpose cri- 
terion. It is also shown that the product states introduced by Wootters in [ p~C| ] form BSA ensemble 
for these states. 

Considering a real matrix parameterization of density matrices, Verstraete et al. in fl4] , |l5| 1 
showed that an local quantum operations and classical communications (LQCC) corresponds to left 
and right multiplication by a Lorentz matrix. They also have shown that the real parameterization 
matrices R{j = Tr(p<Ji<jj) can be decomposed as R = L\ Lj, where L\ and L\ are proper Lorentz 
transformations and nas either a diagonal form correspond to Bell decomposable states or it 
has a non-diagonal form. For the case that the matrix becomes Bell decomposable state we 



presented optimal decomposition for states obtained from them via some LQCC actions [13|. In 



other cases that Yl has non-diagonal form, corresponding matrix is a special cases of ICD states, 
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where we consider it in this paper, and investigation of action of LQCC on these states will appear 
in another paper. 

The paper is organized as follows. In section 2 we introduce ICD states, then we study their 
separability properties. The concurrence of these states is evaluated in section 3, via the method 
presented by Wootters in [1C]. In section 4 we obtain L-S decomposition of these states and we 



prove that thus obtained decomposition is optimal. The paper is ended with a brief conclusion in 
section 5. 



2 Iso-concurrence decomposable states 

In this section we define iso-concurrence decomposable (ICD) states, then we give their separability 
condition. The iso-concurrence states are defined by 

|Vi) = cos0|TT)+sin0|U}), (2-1) 

|V> 2 > = sin0|TT>-co S 0|U>), (2-2) 

Hz) =™s9 HI) + sin0 HI)), (2-3) 

|V> 4 > = sm9 |U)- cos 0|jT»- (2-4) 

It is quite easy to see that the above states are orthogonal and thus they span the Hilbert space of 
2 <S> 2 systems. Also by choosing 6 = f the above states reduce to Bell states. Now we can define 
ICD states as 

4 4 

P = ^Pi\i/>i){il>i\> 0< P i<l, 5>i = l. (2-5) 
i=i i=i 

These states form a four simplex (tetrahedral) with its vertices defined by p\ = 1, P2 = 1, Pz = 1 
and p4 = 1, respectively. 
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A necessary condition for separability of composite quantum systems is presented by Peres || . 
He showed that if a state is separable then the matrix obtained from partial transposition must be 
positive. Horodecki family R have shown that Peres criterion provides sufficient condition only for 
separability of mixed quantum states of dimensions 2 (g> 2 and 2§3. This implies that the state 
given in Eq. (|2-5|) is separable if and only if the following inequalities are satisfied 





-P2) 


< < 


y + 




2 sin 2 20, 


(2-6) 


(P2 


-Pi) 


< 1 


y 4p 3 p 4 + (p 3 




2 sin 2 20, 


(2-7) 


(P3 




< 1 


y 4pip 2 + (pi 


~P2^ 


2 sin 2 20, 


(2-8) 



(P4 - Pa) < V 4piP2 + (pi - P2) 2 sin 2 20. (2-9) 



Inequalities (2-6) to (2-9) divide tetrahedral of density matrices to five regions. Central regions, 
defined by the above inequalities, form a deformed octahedral and are separable states. In four 
other regions one of the above inequality will not hold, therefore they represent entangled states. 
In the next sections we consider entangled states corresponding to violation of inequality (|2-6[) i.e. 
the states which satisfy the following inequality 



(Pi - P2) > V 4 P3P4 + (P3 - P4) 2 sin 2 20. (2-10) 
All other ICD states can be obtain via local unitary transformations. 

3 Concurrence 

In this section we first review concurrence of mixed states. From the various measures proposed to 
quantify entanglement, the entanglement of formation has a special position which in fact intends 



to quantify the resources needed to create a given entangled state H. Wootters in [10] has shown 
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that for a 2-qubit system entanglement of formation of a mixed state p can be denned as 

1 1 



E(p) = H (_ + -VT=C*), (3-11) 
where H (x) = —x In x — (1 — x) In (1 — x) is the binary entropy and the concurrence C (p) is defined 

by 

C(p) = max{0, Ai - A 2 - A 3 - A 4 }, (3-12) 
where the Aj are the non-negative eigenvalues, with Ai being the largest one, of the Hermitian 



matrix R = sfpp^fp and 

p = (a y <g> 0-y)p*((Ty <8> Vy), (3-13) 

where p* is the complex conjugate of p when it is expressed in a fixed basis such as{||t),|t|)},{||t),|||)} 
and o y represent Pauli matrix in local basis {||) , ||)} . 

In order to obtain the concurrence of ICD states we the follow method presented by Woot- 



ters in [ 10 1 . Starting from spectral decomposition for ICD states given in Eq. ( |2-5| ) , we define 



subnormalized orthogonal eigenvectors \vi) as 

\vi) = \fp~i\i>i) , {vi I Vj) =Pi5ij. (3-14) 
Now we can define states \x{) by 

4 

= U *j K> . for i = 1, 2, 3, 4, (3-15) 

j 

such that 

(xi | xj) = (UTU T )ij = \i5ij, (3-16) 

where ry = {vi \ Vj) is a symmetric but not necessarily Hermitian matrix. To construct \x{) we 
consider the fact that for any symmetric matrix r one can always find a unitary matrix U in such 
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a way that Aj are real and non-negative, that is, they are the square roots of eigenvalues of tt* 
which are same as eigenvalues of R. Moreover one can always find U such that Ai being the largest 



one. 



By using the above protocol we get the following expression for state of p given in Eq. ( |2-5| ) 



—pi sin 29 \/P\P2 cos 29 
^fP\P2 cos 29 p2 sin 29 


Now it is easy to evaluate Aj which yields 



V 







P3 sin 29 — y/psPi cos 29 

- y/psPi cos 29 —p A sm29 



(3-17) 



Ai 



\ ((pi ~ P2) sin 29 + ^Ap l p 2 + {p 1 -p 2 ) 2 s\n 2 29^j , 
A 2 = 5 ( (P2 - Pi) sin 29 + J 4p lP2 + (pi - p 2 ) 2 sin 2 29) , 



A3 = 5 ((pa ~ Pa) sin 29 + yj 4p 3 p A + (p 3 - p A ) 2 sin 2 26^ , 



(3-18) 



A 4 = \ ( (p 4 - pa) sin 29 + J Ap z p A + (p 3 - p 4 ) 2 sin 2 26* 



Thus one can evaluate the concurrence of ICD states as 



C = (pi — p 2 ) sin 29 - J 4p 3 p A + (p 3 - p 4 ) 2 sin 2 20. 



(3-19) 



It is worth to note that thus obtained concurrence is equal to the amount of violation of inequality 



(12-101) . 



Finally we introduce the unitary matrix U which is going to be used later 



U 



ict\ 


— ict 2 










a 1 














"3 


— «4 








ia A 


ia 3 



(3-20) 
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((pi + p 2 ) sin 20 + \J 4pip 2 + (pi - P2) 2 sin 2 20 



CH2 



V2 ^4p lP2 cos 2 29 +{pi+p 2 ) 2 sin 2 29 + (p a + p 2 ) sin 20 ^ 4p x p 2 + (pi - P2) 2 sin 2 29 

^J2p\p2 cos 20 



"3 



4pip2 cos 2 20 + (pi + P2) 2 sin 2 20 + (pi + p 2 ) sin 20 J 4pip 2 + (pi — P2) 2 sin 2 20 



((P3 + P4) sin 20 + sj ^P3P4 + (P3 - P4) 2 sin 2 20 



(3-21) 



\/2 4p 3 p 4 cos 2 20 + (p 3 + p 4 ) 2 sin 2 20 + (p 3 + p 4 ) sin 20 J 4p 3 p 4 + (p 3 - p 4 ) 2 sin 2 20 



^/2p 3 p 4 cos 20 

a 4 — 



4p 3 p 4 cos 2 20 + (p 3 + p 4 ) 2 sin 2 20 + (p 3 + p 4 ) sin 20^/4p 3 p 4 + (p 3 - p 4 ) 2 sin 2 20 

4 Lewenstein-Sanpera decomposition 

According to Lewenstein-Sanpera decomposition [JO]], any 2-qubit density matrix p can be written 



as 



P = ^Pse P + (1 — A) \if}) , A G [0, 1] , (4-22) 

where p sep is a separable density matrix and is a pure entangled state. The Lewenstein-Sanpera 
decomposition of a given density matrix p is not unique and, in general, there is a continuum set 
of L-S decomposition to choose from. The optimal decomposition is, however, unique for which A 
is maximal and 

p = + (1 - A(°^)|V (opt) }( V (opt) l , A (op4) G [0, 1]. (4-23) 



Lewenstein and Sanpera in [11| have shown that any other decomposition of the form p = Xp sep + 



(1 - X)\tp)(tjj\ with p / p( op ') necessarily implies that A < A( op ') Q. We refer to Eq. (|4-23|) as 
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optimal decomposition in the sense that A is maximal and p s is the best separable approximation 
(BSA). 

Here in this section we would like to obtain L-S decomposition for ICD states. Let us consider 
entangled state p belonging to entangled region defined by Eq. (2-10). We start by writing density 



matrix p as a convex sum of pure state \tp\) and separable state p s as 

p = Xp s + (1 — A) \ipx) (ipi\ . (4-24) 

where, p s corresponds to intersection of the line passing through the points corresponding to 
and p with 2-dimensional manifold defined by the equality in relation ( |2-6| ) . Expand- 
ing separable state p s as p s = J2i=iPi IV'i) (V^l an d using Eq. fl2-5|) for p we arrive at the following 
results 

p[ = Pl + *~\ P'i = j for i = 2,3,4, (4-25) 

where 



A = 1 - pi + P 2 + J + (pa - PA?- (4-26) 
V sin 29 



Now using the fact that for pure states defined in Eqs. ( |2-1| ) to (2-4) we have C(\tpi)) = sin 29, one 



can easily evaluate average concurrence of the decomposition given in Eq. ( 4-24| ) as 



1 - A)C(|^i)) = (pi - pa) sin 29 - J4p 3 p 4 + (p 3 - p 4 ) 2 sin 2 29. (4-27) 



Equation ( |4-27| ) shows that average concurrence of p in the L-S decomposition given in (4-24) is 
equal to its concurrence. As the concurrence of a mixed state is defined as the minimum of the 
average concurrence over all decompositions of the state in terms of pure states [|10|, so according 



to Eq. ( |4-24| ) we should have C(p) < (1 - X)C(ip). But Eq. Q4-27D shows that the optimal 



decomposition of ICD states saturate this inequality. 
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In the rest of this section we will prove that the decomposition ( 4-24j ) is the optimal one. To 
do so we have to find a decomposition for p s in terms of product states |e a , f a ), i.e. 



Ps — ^ ' A Q \ e a , f a ) (e a , f a 



(4-28) 



such that the following conditions are satisfied [11] 

i) All A a are maximal with respect to p a = A a \e a , f a ) (e a , f a \ + (1 — A) {ipi] and projector 

Pa — | &a j fa) i fa \ • 

ii) All pairs (A Q , A^) are maximal with respect to p aj3 = A a \e a , f a ) {e a , f a \+Ap \ep, fp) {ep, fp\ + 
(1 — A) and the pairs of projector (P a , Pp). 



Then according to 11] p s is BSA and the decomposition given in Eq. ( 4-24 ) is optimal. 



p 1 \ip) 1 > 0. They have also 



Lewenstein and Sanpera in [11] have shown that A is maximal with respect to p and P 
iff a) if \tp) K{p) then A = 0, and b) if \ip) G K(p) then A = 
shown that a pair (Ai, A2) is maximal with respect to p and a pair of projectors (Pi, P2) iff: a) if 
^2) do not belong to TZ(p) then Ai = A2 = 0; b) if \ipi) does not belong, while \ip2) G T^(p) then 
A x = 0, A 2 = (V2I p- 1 c) if IV2) G K(p) and p~ x IV2) = then A; = (^| p~ l \^)~\ 

i = 1, 2; d) finally, if , \ifa) G ft(p) and p~ l |^ 2 ) / then 



(4-29) 



Ai= ((V2I p- 1 1^2) - I (Vilp- 1 !^) |)/A 
a 2 = ((H p' 1 - \ (Hp' 1 ^) \)/d, 

where D = p~ l (^2) P _1 ^2) - I (ipl\p-l [^2) | 2 - 

Now we return to show that the decomposition given in ( 4-24| ) is optimal. Wootters in [1C] has 
shown that any 2 ® 2 separable density matrix can be expanded in terms of following product states 



21 



^ \x x ) + e 102 \x 2 ) + e ie * \x 3 ) + e ie * \x£ 



(4-30) 
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Z2 



*3 



\ z i) 



e^ 1 |xi) + e iea |x 2 ) - e id3 \x 3 ) - e i6i \x 4 ] 



e iei \xi) - e i02 \x 2 ) + e i6z \x 3 ) - e i04 \x 4 ] 



e iei \xi) - e i02 \x 2 ) - e i9z |x 3 ) + e i9i \x 4 ] 



(4-31) 
(4-32) 
(4-33) 



provided that Ai — A2 — A3 — A4 < 0. Now, the zero concurrence is guaranteed by choosing phases 
8i, i = 1, 2, 3, 4 satisfying the relation Xw=i e 2l6j ^j = 0. 

Now using the fact that for marginal separable states p s (located at the boundary of separable 
region) the eigenvalues Aj satisfy the constraint Ai — A2 — A3 — A4 = 0, we can choose the phase 
factors 6i as 9 2 = #3 = #4 = &i + § • Choosing B\ = we arrive at the following product ensemble 
for p s 

\z\) = +a4)\/pilV'i) - a4)\/p£l^2> 

-(a' 4 + ia' 3 )^\i) 3 } - (a 3 - ia 4 ) ^Jp\ |V>4>), 
|z 2 ) = |(— + <4)'\/K \ipi) -i(a[ - oQyf^\ip 2 ) 

V V (4-34) 

N) = |(— «(o^ - a 2)\fp\ +*(o / i + a / 2)\/P2"lV'2) 

+(q' 4 - ia / 3 )y / ^|V'3) + (a' 3 + ia' 4 )^\ip 4 }), 

\z 4 ) = -«2)y5iW +«W + a'2)\fiW 

-(a 4 - ia' 3 )^ z ^3) - («3 + m 4 )-y / ^|'i/?4)), 
where are defined in the same form as given in Eq. ( |3-21 ) provided that they are expressed in 

terms of p[ (coordinates of p s ). 

Now in order to prove that the decomposition given in Eq. ( ggp is optimal we first consider the 

cases that p has full rank. Let us consider the set of four product vectors {\zi)} and one entangled 



state \tpi)- I n Ref. |1C] it is shown that the ensemble are linearly independent, also it is 
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straightforward to see that three vectors \z a ) , \zq) and \tp\) are linearly independent. Now let us 
consider matrices p a = A a \z a ) (z a \ + (1 — A) |-0i) (V>i|- Due to independence of \z a ) and l^i) we can 
deduce that the range of p a is two dimensional, thus after restriction to its range and defining dual 



basis \z a ) and 



■01 ) we can expand restricted inverse p a 1 as p a 1 = A a 1 1 z a ) ( z a | + ( 1 — A) 1 1 V 7 ! ) ( V 7 ! | 



(see appendix). Using Eq. ( |5-35 ) it is easy to see that ^z^l^^^l^o^ — A a ^. This shows that A a are 



maximal with respect to p a and the projector P a = \z a ) (z a \. 

Similarly considering matrices p a p = A a \ z a ) (z a \+Ap \ zp) {zp\ + (l— A) |0i) (-01 1 and considering 
the independence of vectors \z a ), \zp) and we see that rang of p a p is three dimensional where 
after restriction to its range and defining their dual basis \z a ), \zp) and ^ \ we can write restricted 
inverse p~p as p~p = A~ 1 \z a )(z a \ + A^ 1 ^ )( zp\ + (1 — A)~ 1 |0i )(0i|. Then it is straightforward 
to get (e a \ p~p \z a ) = A" 1 , (zp\ p~^ \zp) = A" 1 and (z a \ \zp) = 0. This implies that the pair 
(A a ,Ap) are maximal with respect to p a p and the pair of projectors (P a ,Pp), thus complete the 
proof that the decomposition given in Eq. ( pT22D is optimal. 



We now consider the cases that p has not full rank. First let us consider p' 2 = case. In this case 
the pairs {\z\) , l^)} and also {\zs} , {z^}} are no longer independent with respect to IV'i)- in this 
cases we can easily evaluate |?/>i) in terms of \zi) and \z2) then matrix p\2 can be written in terms 
of two basis \z\) and \z%) which yields after some calculations, (•ZilPi^l 2 -!^ ~~ A2 ("i +a 2) Pi+i 1 A ) 



r i2 

**\p£\z*) = Al K+^)V 1+ d-A) and ^| p -i|^ = ^_A) where Fi2 = a 1 A 2 (o / 1 + a' 2 )Vi + (1 - 
A)(Ai + A2). Using the above results together with Eqs. ( [4-29 ) we obtain the maximality of pair 



(Ai, A2) with respect to p\2 and the pair of projectors {Pi,P%). Similarly we can express in 
terms of 1 2^3) and \z/£) and evaluate p%l, which get (zs\p&t\zs) = A4 ( a i "2) Pi+i 1 A ) ^ /zi\p^l\z4 



A 3 K-^)V 1+ (i-A) and ^ 3 | p -i|^ 4 ^ = dg± wi th r 34 = A 3 A 4 K - a' 2 )Vi + (1 - A)(A 3 + A 4 ), 
together with Eqs. ( [4-29 ) we obtain the maximality of pair (A3, A4) with respect to ^34 and the 
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pair of projectors [P^^PjC)- For other choices of a and j3 three vectors \z a ) ,\zp) and l^i) remain 
linearly independent, thus with the above mentioned method we can prove the maximality of pairs 
(A a ,A^). 

Also in the cases of p' 3 = or p' 4 = 0, three vectors \z a ) ,\zp) and |V>i) remain linearly independent 
thus maximality of pairs (A a , A^) can be proven in the same method of full rank one. 

Now let us consider the cases that rank p is 2. In this case we have to write the separable 
part in terms of rank one projection operators which are not product any more but as far as the 
maximality of pairs (the requirement for the maximality of separable part of density matrix in this 
case) are concerned the separability of the projectors does not play any role. 

For p2 = P3 = we have p' x = p' A = ^, therefore p s can be written as 

Ps = \ (l^i) + |V>4) <V>4|) = \ (|ei, /i> (ei, h\ + |e 2) h) (e 2 , h\) , 

where |ei s 2, fi,2) = (cos 9 |f ) T i sin 9 \ \)) <g> -j= ± i |J,)). Now by writing in terms of product 
states \e%, f\) and |e2, / 2 ) the matrix pi 2 can be expressed in terms of two basis \e\, f±), |e 2 , / 2 ) and 
we get (eij/ilp^kij/i) = 2A2 +^~ A) , ^e 2 , / 2 |pr 2 1 l e 2, ^) = ^rla"^ and ( e ii h\Pii\ e 2, h 

-(l-A) 



where Ti 2 = 2AiA 2 + (1 — A)(Ai + A 2 ). Using the above results together with Eqs. (4-2 



we obtain the maximality of pair (Ai, A 2 ) with respect to p\i and the pair of projectors (Pi,P 2 ). 
The other rank-2 cases can be proved via similar method. 

5 Conclusion 

We have derived Lewenstein-Sanpera decomposition for iso-concurrence decomposable states. It is 
shown that for these states the average concurrence of L-S decomposition is equal to the concurrence 



of the states. It is also shown that Wootters product states defined in [1C] form the best separable 
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approximation ensemble for these states. 
Appendix 

Let us consider the set of linearly independent vectors {|<fo)}, then one can define their dual 
vectors { 4>ij} such that the following relation 

k \<t>jj = Sij (5-35) 



hold. It is straightforward to show that the and their dual { (pi)} posses the following 

completeness relation 

J2\fr)(c/>i\ = i, £i<m<&i = j- ( 5 - 36 ) 

i i 

Consider an invertible operator M which is expanded in terms of states \<pi) as 

M = ^o ij \<l> i )(<f> j \ (5-37) 

i 

Then the inverse of M denoted by M _1 can be expanded in terms of dual bases as 

M' 1 = Y,bSi)(h\ (5-38) 

i 

where bij = (A^ 1 )^ and Aij = aij. 
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